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Abstract. We consider the 2m-th order elliptic boundary value problem Lu = f{x,u) on 
a bounded smooth domain C with Dirichlct boundary conditions u = ■§;;u ~ ... = 
(^)'"~^it = on dfl. The operator L is a uniformly elliptic operator of order 2m given by 
1/ = ( - '^iji^) dx^xj )"" + S|Q|<2m-i baix)D°'. FoT the nonlincarity we assume that 

linis^oo "^^^j^ = h{x), lim^-t-oo ^j^fi^ — k{x) where h,k G C(r2) are positive functions and 
g > 1 if < 2m, 1 < q < j^j;;^™ if A^ > 2m. We prove a priori bounds, i.e, we show that 
||m||^oo(-q-| < C for every solution u, where C > is a constant. The solutions are allowed to 
be sign-changing. The proof is done by a blow-up argument which relies on the following new 
Liouville-type theorem on a half-space: if m is a classical, bounded, non-negative solution 
of {—A)™u = in M.^ with Dirichlet boundary conditions on dR'!^ and q > 1 if N < 2m, 
1 < 9 < #35^ if A^ > 2m then u = 0. 



1. Introduction 

A priori bounds for solutions of elliptic boundary value problems have been of major 
importance at least as far back as Schauder's work in the 1930s. In this paper we prove a 
priori estimates on bounded, smooth domains Q C for solutions of higher order boundary 
value problems of the form 



(1.1) Lu = fix, u) vhVL, u = -— u = . . . = -— n = on dVt. 

Here i' is the unit exterior normal on dQ and 

i,j=l ■' \a\<2m,-l 

is a uniformly elliptic operator with coefficients G and a^j G C^"*~^(f2) such that 

there exists a constant A > with A~^|^p < YlZ=i^ii^^'^^i^i - ^'^'^ ^ G M^, x G fi. 

Our main result is the following: 



Theorem 1. Suppose f2 C M is a bounded domain with G C . Let m G N and assume 

N+2m 
N-2m 



that q > I if N < 2m and 1 < q < jr^^ if N > 2m. Suppose further that there exist 
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positive, continuous functions k, h : Q {0, oo) such that 

(1.2) lim i^^ = h{x), lim l^^ = k{x) 

uniformly with respect to x & fl. Then there exists a constant C > depending only on the 
data aij, b^, ^, N, q, h, k such that ||m||oo < C for every solution u of (1.1). 

Remark 2. Suppose the nonlinearity depends on a real parameter \, i.e, : ^2 x M — R 
and 

lim '^^ = h{x), lim '^'^.^.^i = k{x) 

s^+oo \s1 s-*-oo A|s|^ 

uniformly with respect to x E Vt and AG [Aq, oo) where Aq > 0. Then the a priori hound of 
Theorem 1 depends additionally on Aq hut not on A. This is important in the study of glohal 
solution branches of a parameter dependent version of (1.1), which we will pursue in future 
work. 

We are focusing on the case of superlinear nonlinearities f{x,u) with subcritical growth. 
A model nonlinearity is f{x, s) = {sl*^. Our results hold with no restriction on the shape of 
the domain Q and for general, possibly sign-changing solutions. This is important since the 
lack of the maximum principle for higher order equations does not allow to restrict attention 
to positive solutions only. 

In the second-order case m = 1 a priori bounds for positive solutions have been established 
for subcritical, superlinear nonlinearities via different methods by Brezis, Turner [7], Gidas, 
Spruck [13], DeFigueiredo, Lions and Nussbaum [10] and recently by Quittner, Souplet [20] 
and McKenna, Reichel [19]. In the higher-order case m > 2 the theory is far less developed 
and strongly depends on the type of boundary conditions considered. For Dirichlet boundary 
conditions we only know of a result of Soranzo [23] , who proved a priori bounds for positive 
radial solutions on a ball if L = (—A)™. For Navier boundary conditions the picture is more 
complete. Let L = (— Lq)™ where Lq = Qij q^q^ + ^a^f" is a second order operator and 
suppose the boundary conditions are of Navier-type: 

(1.3) u = {-L)u = ... = {-1)""-^ = on an. 

Soranzo [23] proved a priori bounds for positive solutions if Lq = A and is a bounded 
smooth convex domain. Recently, Sirakov [22] improved this result to general operators 
L = (— Lq)™ and general bounded smooth domains. Both authors strongly use the fact that 
the boundary conditions (1.3) allow to write the problem as a coupled system of second order 
equations, where each equation is complemented with Dirichlet boundary conditions. In this 
case maximum principles are available. In contrast, the higher order Dirichlet problem can 
not be rewritten as a system and therefore requires different techniques. 

In our approach we extend the so-called "scaling argument" of Gidas and Spruck [13], 
which they used to deal with the second order case m = 1 and positive solutions. Let us 
give a brief sketch of their method. Gidas and Spruck assume that there exists a sequence 
of positive solutions with L°°-norm tending to +oo. After rescaling the solutions to norm 
1 and blowing-up the coordinates one can take a limit of the rescaled solutions and obtains 
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a nontrivial positive solution of a limit boundary value problem —Au = u'^ on either the 
full-space or the half-space = {x G : xi > 0} together with Dirichlet boundary 
conditions. Then a contradiction is reached provided that a Liouville-type result is available, 
i.e., a result which shows that the non- negative solutions of the limit problem must be 
identically zero. For subcritical q Gidas and Spruck [13], [14] proved both the full-space and 
the half-space Liouville theorem for —Am = u'^ via the method of moving planes. 

In order to deal with the higher order Dirichlet problem (1.1) and solutions which may 
change sign, the blow up procedure has to be modified. Indeed, even under assumption 
(1.2), there seems to be no direct argument to exclude the case of negative blow up (i.e., the 
existence of a sequence of solutions which is not uniformly bounded from below). Instead, it 
is excluded a posteriori after passing to the limit equation. Once this is done, we still need 
Liouville theorems for nonnegative solutions of the higher order problems on M^, M.^. The 
full-space Liouville theorem stated next is already known; it was proved by Lin [18] if m = 2 
and for general m > 2 by Wei, Xu [24]. 

Theorem 3 (Wei, Xu). Let m G N and assume that q>lifN<2mandl<q< if 
N > 2m. If u is a classical non-negative solution of 

(-A)"n = n« m M^, 

then M = 0. 

Even in the case of the Navier boundary conditions, the corresponding Liouville theorem 
for the polyharmonic problem in the half-space is harder to achieve and has only recently been 
proved by Sirakov [22]. Due to the lack of a (local) maximum principle, the corresponding 
Dirichlet problem is even more difficult to deal with. Here we show the following new Liouville 
Theorem for the half-space which complements Theorem 3. 

Theorem 4. Let m G N and assume that q > 1 if N <2m and 1 < q < if N > 2m. 

If u is a classical non-negative bounded solution of 

o Cynx — 1 

(1.4) (-A)™M = u'^ m Mil, u = —u = ... = -ru = on dR^ 

oxi ax™ 

then u = 0. 

We point out that the critical case q = ^'tlm allowed in Theorem 4. Note also that 
Theorem 4 holds in the class of bounded solutions. It remains an open problem to extend 
the result to the class of all (possibly unbounded) classical positive solutions. 

Let us outline the proof of Theorem 4 and point out the main difficulties. Following Gidas 
and Spruck, we transform the half-space problem via a Kelvin inversion into a problem in the 
unit ball, where the point at infinity is mapped onto the boundary point P = (—1, 0, . . . , 0). 
The transformed solution satisfies 

{-A)"'v = 22'"|x - P|-"t;'? pointwise in 5i(0), 

m— 1 

v = on 95i(0) \ {P}, 



_d__ _ fd_\ 
dv ' ' ' \di' J 
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where a = N +2m—q{N —2m) > 0. The key step is to show that v is axially symmetric around 
the xi-axis. In the second order case m = 1 this is proved with the classical moving plane 
method, which is a local method based on the maximum principle. The same local approach 
fails for the higher order case m > 2 since the maximum principle is not available. Very 
recently, a new development in the moving plane procedure by Berchio, Gazzola and Weth [4] 
overcame part of this difficulty. The authors made the moving plane method applicable for 
classical solutions of polyharmonic Dirichlet problems on balls. Instead of a local maximum 
principle method they argue via the Green integral-representation and properties of the 
Green function. However, the method of [4] does not apply here, since the solution v of 
(1.5) may have a singularity at P G dBi{0). To overcome this problem, a large part of 
this work is devoted to show that every solution v of (1.5) which corresponds to a bounded 
solution of (1.4) can be of represented via the Green function. In this step we also use Green 
function estimates of Grunau and Sweers [15]. Then we apply a moving plane argument, 
using the Green function representation and the Hardy-Littlewood-Sobolev inequality, to get 
the desired symmetry result. Comparing this variant of the moving plane method with the 
one in [4], we point out that Berchio, Gazzola and Weth allow more general (non-Lipschitz) 
nonlinearities, but their argument relies on Green function representations for directional 
derivatives of the solution which in our situation might not exist. 

Once the symmetry result for v is established, we readily conclude - following Gidas and 
Spruck [13], [14] again - that the corresponding solution u of (1.4) is axially symmetric 
around any axis parallel to the xi-axis. Consequently, m is a function of xi only and hence 
solves an ordinary differential equation. It is then easy to conclude that u = 0. 

We recall that the original moving plane method goes back to Alexandrov [3] and Serrin [21] 
and was further developed by Gidas, Ni, Nirenberg [12] for second order equations. Recent 
improvements of the moving plane method for higher order equations and pseudo differential 
operators using integral representations rather than local maximum principles were achieved 
by Chang, Yang [8], Berchio, Gazzola, Weth [4], Li [16], Chen, Li, Ou [9] and Birkner, 
Lopez-Mimbela, Wakolbinger [5]. 

The paper is organized as follows. In Section 2 we prove Theorem 1 assuming Theorem 4. 
We give the details of the blow-up procedure taking into account that we allow for solutions 
blowing up to either at +oo or — oo. The rest of the paper is devoted to the proof of Theo- 
rem 4. In Section 3 we prove a Green-representation formula on half-spaces (Theorem 9) by 
approximating the half-space by a family of growing balls. Based on the Green-representation 
for balls and by a careful estimate of the boundary integrals and the monotone convergence 
theorem we obtain a Green-representation for the half-space. Finally, in Section 4 we prove 
Theorem 4. 



2. Proof of Theorem 1 - the blow-up argument 

In this section we give the details of the blow-up argument for the proof of Theorem 1 
under the assumption of the validity of the Liouville-type result of Theorem 4. The proof 
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uses standard linear U^-W"^^'^ estimates for linear problems 



(2.1) 



Lu = g{x) in Q, 



(2.2) 




du 



) 



M = on dVL. 



Recall the following basic estimate of Agmon, Douglis, Nirenberg [2]. 

Theorem 5 (Agmon, Douglis, Nirenberg). Let Vt C be a hounded domain with dfl G 
C2"\m e N. Let a.,^ E C^'^-^m, 6„ G L°°(fi), g e LP(fi) for some p E (l,oo). Suppose 
u G W'^'^'^iyt) n W^'^iVt) satisfies (2.1). Then there exists a constant C > depending only 
OTl II ti^j II (72m — 2 ^ ll^alloo, A, r2, A^, p, m and the modulus of continuity of Oij such that 



We will also be using the following local analogue of this result. Though the proof may be 
standard we give it for the reader's convenience. 

Corollary 6. Let n be a ball {x E : \x\ < R} or a half-ball {x eR^ : \x\ < R,Xi> 0}. 
Let m E N, Oij E C2"^-2(n), E L°°{n), g E LP{n) for some p E (l,oo). Suppose 
u E W^'^^PiQ) satisfies (2.1) 

(i) either on the ball 

(ii) or on the half-ball together with the boundary conditions u = -^u = . . . = -^^^^u = 
on {xER^ : \x\ < R,Xi = 0}. 

Then there exists a constant C > depending only on \\aij\\c'2m^2,\\ba\\QQ, N,p,m, the 
modulus of continuity of Oij and R such that for any a E (0, 1) 



Proof. It is sufficient to prove the result for R = 1. For a E (0, 1) let rj E Cq"^{Bi) be 



||w||vK2'".P(Q) < C'(||fl'||LP(Q) + ll^lUp(n))- 



||M||vi/2-.p(f7nB,K) < 



c 



(IblUp(n) + \\u\\lv{u)). 





^^ and 






|/3|<2m-l 
|7|<2m-|/3 
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where cp^^ are bounded functions with ||c/3^^||oo < Ci and Ci = C{m) max{||6Q,||oo, ||ajj||c2™-2}- 
By Theorem 5 

||V^"'u||ip(nnB<,) < C'2(^||5'||lp{Q) + ^ \\^^u\\Lp{nr\B^,){'^ - o'Y^'^ 

0<fc<2m-l 
Q<l<2m-k 

2m- 1 

<C,{\\9\\Lnn) + Yl l|V'n|U.(^nB^,)(l-a)'^-2-). 

k=0 

If we introduce for k E Nq the weighted norm $fc = suPq^^<]^(1 — cr)''|| V'^n||LP(Qns<^) then the 
last inequahty imphes 

2m.- 1 

(2.3) $2m <C3(||(?||LP(f7)+ ^ 

k=0 

Recall the standard interpolation inequality, see Adams, Fournier [1] , for < k < 2m — 1 

II V''n||LP(QnB<,) < e||V^"'u||LP(QnB^) + Ce^^\\u\\LP(nnB„), 

where C is homothety invariant and hence independent of a. Using this we find that for 
every fixed 6 > there exists cr(5) e (0, 1) such that 

< (1 -0^)1|V'=M||LP(nnB.) +^ 

kf /I _N2m— fc , II T-72m„ . II i , „ , / 1 fe 



< (1 - - a)^™-\|| V^'"iz|U.(f,nB.) + Ce5^(l - a)-'=||M|Up(f,nB.) j + 5 

= e(l - fT)2'"||V2™M||LP(nni?.) + Ce^\\u\\LPinnB^) + S. 
Since 6 > was arbitrary, we see that $fc ^ e*^'2m + C'e^o- Hence it follows from (2.3) that 

$2m < C'4(||5'||lp(C) + ||m||lp(!^)), l-G., 

||V^"u||LP(nnB.) < ■^^-r^^(II^IU''(f^) + \\u\\LPin))- 
Using the interpolation inequality again we obtain the claim. □ 
Proof of Theorem 1. It is convenient to rewrite the operator L in the form 

L = {-ir Yl «-(a;)/^"+ Ca(a;)D". 

|a|=2m |a|<2m-l 

Here aa{x) = (^hiii^) ' o.isiii^) ' ' ' (^i2m-ii2mi^) ^ where A^q, is the set of all vectors / = 

(ii, . . . , i2m) G {1, • • • 5 A^}^™ satisfying #{j : = /} = a; for / = 1, . . . , A^. Hence is 
continuous on Q and G L°°{Q). Assume for contradiction that there exists a sequence Uk of 
solutions of (1.1) with := ||Mfc||oo^ooasA;^oo. By considering a suitable subsequence 
we can assume that there exists Xk G Q such that either Mk = Uk{xk) for all /c G N (positive 
blow-up) or Mk = —Uk{xk) for all /c G N (negative blow-up). Define 

1 IziR 

My) ■= xr^fc(^fe""2/ + ^fc)- 
Mk 



A PRIORI BOUNDS FOR HIGHER-ORDER ELLIPTIC PROBLEMS 



7 



Then Hi^fcHoo = 1 and either ^^(0) = 1 for all A; G N (positive blow-up) or Vk{0) = — 1 for all 
/c e N (negative blow-up). We may also assume that — > x G fi. 

Case 1: x E Q. In this case Vk is well-defined on the sequence of balls -Bpfe(O) with pk : = 
M^^™ dist(xfc, dil) — > oo as /c — > oo. Note that 



For y e Bp^ (0) let 



(2.4) (t{y) := a,(Mpy + x^), c^(y) := M^' '^^^--'^c„(Mpy + x^] 
and define the operator 

(2.5) L'^ :=(-!)- 5^ a^(y)D"+ ^'^D'^- 

a|=2m |a|<2m— 1 

The function Vk satisfies 



1-9 



(2.6) L'^Vkiy) = fk{y) in i?p,(0), where My) := _/(M y + x^, M,t;,(i/)). 

By our assumption on the nonlinearity /(x, s) we have that ||/fc||L°°(_Bp^(o)) is bounded in k. 
Note that while the ellipticity constant and the L°°-norm of the coefficients of are the same 
as for L, the modulus of continuity of is smaller than that of Oq,. By applying Corollary 6 
on the ball -B_r(0) for any R> and any p > 1 there exists a constant Cp^R > such that 

\\vk\\w^"--piBii{o)) < Cp^R uniformly in k. 

For large enough p we may extract a subsequence (again denoted Vk) such that Vk — > v 
in C'^'^~^'°'{Bii{0)) as k ^ oo for every R > 0, where v e Cfj^~"^'"(]R^) is bounded with 
IK' I loo = 1 = ±t'(0). Taking yet another subsequence we may assume that fk^F in L°°{K) 
as A; ^ oo for every compact set K C M^. Also we see that 

h{x)v{yY if v{y) > 0, 
kix)\v{y)\'^ ift;(y) <0, 



(2.7) Fiy) 



because, e.g., if v{y) > then there exists ko such that Vk{y) > for A; > /cq and hence 
MkVk{y) — oo as — i> oo. Therefore the assumption on /(x, s) implies that fk{y) — ^ 
h{x)v{yY a.s k —* oo, and a similar pointwise convergence holds at points where v{y) < 0. 
Finally, note that the pointwise convergence of fk on the set = G : v{y) > 0} 
and Z~ = {y G : v{y) < 0} determine due to the dominated convergence theorem the 
weak*-limit F of fk on the set Z+ U Z^. Since b'^{y) — > and a^{y) — > ^^(a;) as /c ^ oo 
and since we may assume that Vk v in W^^^{M.^) we find that v is a bounded, weak 
W^i"/(K'^)-solution of 

-^aij(x)— — -j . 

\a\=2m 1,3=1 ^-^ 
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Since F G L°°(M^) we get that v G iyf„^'^(M^) n is a bounded, strong solution 

of (2.8). Because D'^'^v = a.e. on the set {y G : v{y) = 0} we see that u is a strong 
solution of 

r h{x)v{y)'' ifv{y)>0, 

Cv = I if v{y) = 0, 

[ k{x)\v{y)\'^ iiv{y) < 

in M^. Notice that the right-hand side of the equation is C^(M^). Hence f is a classical 
Cf^'°'{M.^) solution. By a linear change of variables we may assume that v solves 

h{x)s'^ if s > 0, 



(2.9) (-A)™t; = g{v) in M^, where g{s 



k{x)\s\'^ ifs<0. 



By Lemma 15 of the Appendix we find that v > 0. This already excludes negative blow-up 
and implies that g{v{y)) = h{x)v{yY^ t'(O) = 1. Theorem 3 tells us that this is impossible. 
This finishes the contradiction argument in the first case. 

Case 2: x G dVt. By flattening the boundary through a local change of coordinates we may 
assume that near x = the boundary is contained in the hyperplane Xi = 0, and that Xi > 
corresponds to points inside Q. Since dQ is locally a C^'"-manifold, this change of coordinates 
transforms the operator L into a similar operator which satisfies the same hypotheses as L. 
For simplicity we call the transformed variables x and the transformed operator L. Now the 

9-1 

function Vk is well-defined on the set -8^^.(0) Ci {y E : yi > —Mf^^"" Xk^i}- Since 
1 = I VkjO) - Vk{-M^Xk,u 0, • • • , 0)J < M^Xk^VvkWoo 

=±1 =0 

£-1 

we see that either M^^™ Xk^i is unbounded and we can conclude as in Case 1, or (by extracting 

a subsequence) Tk ■= M^"^ Xk,i ^ r > as A; ^ oo. In this case we make a further change of 
coordinates and define 







■= Vk( 




- Tk-, 


^2,- 


■ ■ yZN 








(^1 


- Tfc 


Z2i ■ 


■ ■,zn 








>1 


- Tfc, 


Z2, ■ 


..,Zn 



and likewise the operator L^. Note that Wkirk, 0, . . . , 0) = ±1. Let = G : 2i > 0} 
and = -Br(O) fl for i? > 0. For k sufficiently large the coefficients a^, and the 
operator L'' are well-defined in B^. As before Wk satisfies 

L^Wk{z) = fk{z) in 5^, where fk{z) := ^/(M^"'" z + (0, Xk,2, • • • , Xk,n), MkWk{z)). 

together with Dirichlet-boundary conditions on {z G : \z\ < R,zi = 0}. Hence we may 
apply Corollary 6 on the half-ball B^ for any R> and find that for any p > 1 there exists 
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a constant Cp^R > such that 

\\wk\\w2m.p{B+) - ^p,R uniformly in k. 

As in Case 1 we can extract convergent subsequences Wk ^ w in Cfj^~"'^'°(M+) and fk F 
in L'^iM.^) as k ^ oo, where F > 0, ^ is determined in the same way as in Case 1. This 
time, w is a bounded, strong ^^^^"'^(M^) n Cf„';!"^'"(M^')-solution of 

£u; = F in R^, —w = ... = = on dR^ 

ozi dz^ 

with C as in (2.8). By a hnear change of variables we may assume that w solves 

(2.10) (-A)'"^ = g{w) in M^^, —nj = ... = ^-^w = on M^^, 

where g is defined as in (2.9) of Case 1. The representation formula of Theorem 9 shows that 
w is positive and that g{w{z)) = h{x)w{zY. Therefore w is a positive, bounded and classical 
solution C^^-solution of {—^)^w = h{x)w'^ in with Dirichlet boundary conditions and 
w{0) = 1. A contradiction is reached by Theorem 4. □ 

Proof of Remark 2. Take sequences of solutions {uk,Xk) such that Mk := ||tifc||oo ^ oo as 
/c — >• cx), Afc > Ao > and define the rescaled functions 

My) ■= TT'^ki — y + Xk 



Mk \\ Xk 

Due to the assumption > Aq > one has that M^~'^ /Xk ^ as A; ^ oo. Define further 
al{y) ■.= aa\[--^j y + Xkj, 

^fe):=Mr"'S-"Ai^-'c„((^)"% + .. 

with the corresponding operator L^. Then Vk satisfies 

1 / / M"^~'^\ 1/2™ \ 
L'^Vkiy) = fk{y) where fk{y) := [y^J ^ + ^'^^'^•(^)J ' 

Note that limfc^oo fk{y) = h{x)v{y)'^ on Z+ and similarly on Z~ . The rest of the proof is as 
before. □ 

3. Green representation 

The main result of this section is Theorem 9. There we state conditions on a function u 
on the half-space M.^ under which the Green representation formula 



uix] 



Gt,{x,y){-Aru{y) dy hi all x G 
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holds. Here is the half-space Green function, see (3.1) below. In the next section, in the 
proof of Theorem 4, this representation formula will be applied to solutions of (1.4). 

Let us fix some notation. We recall Boggio's celebrated formula [6] for the Green function 
of the operator (—A)'" with Dirichlet boundary conditions on the unit ball B = {a; G : 

|x| < 1}: 

Gi(x,2/) = A;-|x-yp™-^ / ^ JT^dz 

Ji ^ 

for x,y G B. Here is a suitable normalization constant. By dilation we find the Green 
function for the ball Br = {x G : |x| < R} as follows 



(1-1 


x\ 


m 




\y\ 






\x - y\ 


2 



Gn{x,y)=R'-~''G, (|,|) 



{R'-\ 


\x 


n(i?2 




\y\ 






\x - y\ 


2 



Next we set Pr := (i?, 0, . . . , 0) G and we denote by := {x G : \x - Pr\ < R} the 
ball of radius R shifted by Pr. If we let G~^ denote the Green function on B^ with respect 
to Dirichlet boundary conditions then we find the explicit formula 

Uj,[x,y)-K ^il^ ^ ' ^ J - 2 '"^ ^' J, {z + l)^/^"^' 

with 

; + , , {R' - \x - Pr\'){R' - \y ~ Pr\') 

Finally, if we let G^ denote the Green function of the operator (—A)"* on the half-space 
subject to Dirichlet boundary conditions then 

(3.1) G'+(x,y) = %|x-yp-^ / T^-^dz with ^^(x,y) ^"^^^^ 



2 ' ^' Jo {z + 1)^/^ ^ooy-^^yj \x-y\^ 

for X, y G M^. 

Lemma 7. The Green function on B^ converges pointwise and monotonically to the 
Green function on M.^. 

Proof. The pointwise convergence is easily checked. Let x,y E B^. The monotonicity of 
G~^{x,y) with respect to R is equivalent to the monotonicity of ipR{x,y) with respect to R. 
Thus 



d {R^-\x-PR\^){R^-\y-Pj 



R\ 



dR i?2 

^ ^ r>2 l„ r> \2\/r>2 i . r, |2\ , ^Xi 2 1 , o |2\ , / d2 I™ r> l2\ 



^^{R' - |x - Pr^' -\y- Pr\') + ^(i?^ - |y - + ^(i?^ - |x - P. 
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Setting a := {x — Pb)/Ri b := {y — Pjij/R we have |ap, |6p < 1 and we obtain from the 
previous computation 

d (i?2_|^_p^|2)(^2_|^_p^|2) 



dR i?2 
= 2R 



(1 - |ap)(l - \b\') + (1 + ai)(l - m + (1 + - lap)) 

= R{{ai + 1 + \\am - m + + l + Imil - la^)), 

^ v ' ^ V ' 

i(|a|2+2ai+l) i(|fc|2+26i+l) 

and clearly |ap + 2ai + 1 = |(ai + 1, 02, . . . , a7v)P > 0. This establishes the proof. □ 

In [15], Lemma 3.4, Grunau and Sweers proved the following estimates for the polyharmonic 
Green function Gi on the unit ball if |/c| > m and x G B, ?/ G (9B: 

(3.2) \D^yG^{x, y)\ < Ck,N,m\x - y\"--^-\^\{l - \x\r 

for some constant Ck,N,m > 0. For the Green function Gr on Bfi and on the estimate 

(3.2) transforms as follows: 

(3.3) \D';GR{x,y)\<Gk,N,m\x-yr-''-\'KR-\x\r 
if X G Br, y G dB^. Likewise, 

(3.4) \D'yG+ix,y)\<Gk,N,rn\x-yr-''-\'\\xr 
ifx= {xi, 0,...,0) e B+ with xi e{0,R),ye dB+. 

Lemma 8. Let G be the Green function of (—A)™ with Dirichlet boundary condition on an 
arbitrary ball i? C M" with exterior unit normal v on dB. For any function v G C^™~^(i?) fl 
W'^'^'^{B) with p > ^ one has the following Poisson-Green representation for x G B: form 
even 

m/2 

(3.5) v{x) =J2f U'~'v{y)d.A7"'Gix, y) - A--G(x, y)d.^A'-\{y)) dsy 
,=1 JdB ^ ^ 



Gix,y){-Arviy)dy. 



B 



and for m odd 



(m-l)/2 



(3.6) v{x)=- J2 f {A'~My)duAy-'G{x,y)-A^-^G{x,y)d.^A^-'v{y))ds, 

i=l -JdB 

A("-i)/2i;(y)9.,A("^-i)/2G(x,2/)rf.,+ / G{x,y){-Arv{y)dy. 



dB J B 
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Proof. First assume v G C^™(5). Consider the identity 

m m 

^div (A*"^wVA™-^G - A"-*GVA*-^i;) = ^ (A^"^wA™"^+^G - A""*GA^w) 

1=1 i=l 

= t^A^G - GA^'v in S. 

If we integrate this identity over B and take into account that DyG{x, y) = for \a\ < m — l 
and X E B,y ^ dB then we obtain the claim. For v G C'^^~^{B) fl W'^"^'^{B) we can argue by 
approximation and Lebesgue's dominated convergence theorem if we take into account that 
J^G{x,y)\h{y)\dy < const. ||/i||Lp(B) provided h G Lp{B) and p > □ 

Theorem 9. Suppose that u G C2™"1(mJ) n W^^^'^iR^), p > £^ ts a function with the 
following properties: 

(i) u and all partial derivatives of u of order less than or equal to 2m — 1 are bounded, 

(ii) u satisfies Dirichlet boundary conditions on dM.^ , 

(iii) (— A)™M G L^^^{R^) is non-negative in M^. 

Then 



(3.7) u{x) = / G^{x,y){—A)"^u{y) dy for every x E 



Proof. Let us first consider the case where m is even. It clearly suffices to prove (3.7) for 
X = (xi,0, ... ,0) G with xi > fixed. In the following we consider R > 1x\. Then 
X G B^ , X\ G (0, -R), and (3.3) yields for y G dB^ and « < y the following estimates: 



Combining this with (3.5), we get 



Gl{x,y){-ATu{y) dy - u{x) 

m/2 



\x — 


y\ 


—m—N+2i 1 ™ |m 


\x — 


y\ 


-m-N+2i-l\ \m 



<G\xrj2{ (\^"My)\\^ - i/r'"-^+2*-^ + \x- 2/|-'"-^+2*|9.,A-i«(y)i) dsy. 

i=l -^93+ 



Since |x — y| > |x| for y G dB]^, we conclude that 

m/2 



/ GUx,y){-Ay"^u{y)dy-u{x) < C.^^ / |x - y|-^(|A-Vy)| + |9.,A-Vy)|) rfs, 

m/2-1 

(3.8) =c,J2f ^\x-y\-''{\^My)\ + \duAMy)\)dsy. 



i=0 -^931 
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We claim that, for every x G M^, 



(3.9) 



m/2-l 

E 



\x 



i=0 



y\~'' \A'u{y)\ + \d,^A'u{y)\) dsy ^ as R oo. 



For = 1 this is obvious since 



\x - |/r^(|A^'M(y)| + \d,^A'uiy)\) dsy = \x- 2R\~'[\u^^'\2R)\ + \u^'^''-'\2R)\) 



for z < Y — 1 as a consequence of the boundary conditions. For N > 2 and < a < 6 < 2R 
let us consider the set B{a,b) := {y G dB^ : a <yi <b}. On B{a,b) we have y = {yi,y') 
with \y'\ = \j2Ryi — yf. For > 3 we parameterize B{a, b) by the map 

(a,6)x§^-2 ^ 5(a,6), 

(2/1, ^ (2/1, V2i??/i -|/2^((^)j, 

where = ((/?2, • • • , V^tv-i) and 



F : 



/ cos(/92 sincyJs sin(y94 . . . sin(y9Ar_i \ 
sin sin Lp^ sin (^94... sin ^Pn-i 
cos 933 sin 994 . . . sin (y9jv_i 



V 



COS Sm V9Ar_i 

COSV^AT-l / 



Let DF = (61I62I ■ ■ ■ \bN~i) be the Jacobian Matrix of the map F and gr(DF) = det(DF^ ■ 
DF) = det(6i ■ bj)ij=i n^i be the Gram determinant of DF. Since bi = (l, /^"^^ d) 

and 6, = (0, v/2i?yi-|/?^)^ for ^ = 2, . . . , iV - 1 we find 61 ■ 61 = ,^^5;^, 61 ■ 6^ = for 
i = 2, . . . , — 1. Therefore 



v/gr(DF) = v/67^-(2%i 



Since 



we obtain finally 



2/1) ^ 



det(-D6') = (—1)^ ^ sin (y93(sin 794)^ . . . (sin(y9Ar_i 



det(Z^^)| = R{2Ry^ - y\)^ ■ \ det(Z^^)|. 



>Af-3 



Vi^CDF) < R{2Ry^ - y^; 



Therefore we can write the surface integral as follows (if = 2 the line integral is parame- 
terized by (2/1, ±^/2Ryi - yf 



2\T\ 



\x — y\ ^ dsy = < 



b 



\X 



F{y,,ip)\~''^giDFdipdy^ if AT > 3, 



B{a,b) 



2A-1/2 



R{2Ryi - yj 
„ xl + 2Ryi - 2xiyi 



dyi 



UN = 2. 
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Since \x — F{yi, = xf + 2Ryi — 2xi?/i and h < 2R, we can now estimate as follows: 

iV-3 



Bia,t) Ja {xl + 2Ry,-2x,y,)f Ja (| + 2t(l-f))f 

„9 /„ N~:i „n JV--3 ,„ , iV-3 



with C2 = 2^~ci. Here we have also used that R > 2xi and N > 2. From now on we assume 
a < R and split the remaining integral as follows: 

(3.10) \x-y\~''dsy<^ dt + ^ dt . 

" V ' ^ V ' 

= .h =.I2 



B{a,b) 



In the first integral Ii we have (2 — t) > 1 and therefore 



C2 



1 



2 



If a > 0, we conclude that 

(3.11) /i<| !\-'Ut< 



R Jj^ yoR 
while for a = the substitution z = yields 



"1 



(3.12) r-^dt = ^ dz<^ r dz = ^. 



For lo we have 



(3.13) ^^^|/(2-^)"''^^ = l- 

Collecting the inequalities (3.10), (3.11), (3.12) and recalling that R > 2xi, we obtain 

1/xi for a = 0, 



(3-14) f \x-y\-''dsy<c, . , / ^ 

JB{a,b) \l/VaR for a > 0. 

Now let £ > 0. By the standard mean- value theorem using the Dirichlet boundary conditions 
and the boundedness of the derivatives of orders up to m, there exists 6 > such that 
m/2-1 

J2 {\^Hy)\ + \duAMy)\) <e on e : \y,\ < 6}. 

i=0 
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Hence we apply (3.14) and obtain 

m/2-l 



J2 / \x - yr"" (\A"u{y)\ + \d,^A^u{y)\) dsy 

^=0 -^dBl V / 



< e (f) \x — y\ ^ dsy + cq (b \x — y\ ^ dsy 

B{0,S) Jb{5,2R) 
C5 C5C6 C5 

< e 1 — 1= — * e — as it 00. 



xi V6R xi 
Since e was chosen arbitrarily, we conclude that (3.9) holds. 

Using (3.8), (3.9) and Lemma 7 together with the monotone convergence theorem we get 

u{x)=\im[ GUx,y){-Aru{y)dy= [ Gt,{x,y){-Aru{y)dy. 

This finishes the proof in the case where m is an even integer. In the case where m is odd only 
minor modifications are needed. We use (3.6) instead of (3.5) together with the estimates 
arising from (3.4): for x G y G dB^ and i < 

|A™-G+(x,2/)| < Q,^,^\x - 

and for t < 

By essentially the same estimates as before, we now obtain 



GUx,y){-Aru{y) dy-u{x) 



Bt 



(m-l)/2 (m-3)/2 



i=0 i=0 

Again, as a consequence of the boundary conditions, for every e > there exists 6 > such 
that 

{■m-l)/2 (m-3)/2 

J2 |AV(y)l+ E \duyAHy)\<e on {yeR't : \y^\<6}. 

i=0 i=0 

We therefore may conclude as in the case where m is even that 

{m-l)/2 (m.-3)/2 



''^^R i=0 i=0 

Using again the monotone convergence theorem, we conclude that 



u{x)=Yim GUx,y){-Aru{y)dy= G+ (x, y)(-A)-n(y) dy. 
Jb+ Jr^ 



□ 
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4. Proof of the Liouville Theorem in the half-space 
This section is devoted to the proof of Theorem 4. Let m G N and assume that g > 1 if 

N+2m 
N~2m 



AT < 2m and 1 < g < ii N > 2m. Let M be a classical non-neerative bounded solution 



of 



-A)™?/ = in M^, u= — = ... = r = 0on9... 



dxi "' dxl 
We need to show that n = 0. From Theorem 9 we know that 

(4.1) u{x) = I G~^{x,y)u^{y) dy for every x G M:!^. 
We consider the conformal diffeomorphism 

\y + eiY 

where ei = (1, 0, . . . , 0) is the first coordinate vector. The following formula shows how Gji 
is related to the Green function Gi on the unit ball. 

/ \ 2m-N 

Lemma 10. G+((/?(x), (/?(?/))= (^^^r^^^^^j Gi{x,y) foraUx,yeB. 
Proof. An easy calculation yields 

(4.2) |(^(x)-y(y)|= , ,^'"^,7^1 : fora;,yGB. 



Considering the functions ip{x,y) = ^^i-^^p ^^"^ ''Pooi.x.y) = as in Section 3, we 

obtain 

i)^{i^{x),'^{y)) 



\x + ei 


\^ipi{x)\y + ei 


|Vi(2/) _ (2x1 + 2- 


|a; + ei|2)(2yi + 2- 


\y + ei 






\x - y\ 


2 






\x 


- y\ 


2 



(1-1 


x\ 






\y\ 


?) 




\x-y\ 


2 



ip{x, y) for X, y G B. 



We conclude that 



G'i(^(a:),^(y)) = -fW{x)-^{y)\'-~^ (TTT)^'^^ 

2m-Ar ri'{x,y) ^m-l 



a: + ei||i/ + ei|/ 2 ' ^' Jo (;2 + 1)^/' 



[- 1- r) Gi(x, y) for all X, ?/ G B. 

V X + ei -y + ei / 



F + eilly + ei 

Corollary 11. Define the function f : B — > M 6?/ 

w(x) := |x + eip'"-^M(v5(x)) 



□ 
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and the function h : B x [0, oo) —>■ [0, oo) by 

h{x,t) := 2^'^\x + ei\-'^ 
where a := N + 2m — q{N — 2m) > by assumption on q. Then v satisfies 

(4.3) v{x) = / Gi{x,y)h{y,v{y)) dy forallx^'B. 

Jb 

Proof. The Jacobian determinant of (/? satisfies \Jip{y)\ = i^^^^piv for y £ Therefore we 
have 



u{^{x))= I Gt,{ip{x),y)u\y)dy = j^GMx),^{y))u\^{y))\J^{y)\ dy 

^ II ^ I G,{x,y)(\y + er-'"'v{y)) dy 

x + ei\\y + ei\J V / |y + ei|"'^^ 

Gi{x,y)h{y,v{y))dy 



\2m-N 



\X + Ci 

for all X G B, so that 

v{x) = \x + ei\^"'-''u{^{x)) = I Gi{x,y)h{y,v{y))dy 



□ 

Proposition 12. The function : B ^ M zs axially symmetric with respect to the Xi-axis. 

Proof. We assume that N > 2 and v ^ 0, since otherwise the statement is trivial. The 
integral representation (4.3) implies that v is strictly positive in B. Note that for every 
Xq G 9B \ {— ei} we have that 

lim v{x) = lim \x + ei |2'"-^u((^(x)) = |xo + ei ^'"-^u f 2^^^^ - ei) = 
x^xo,xeB x^xo,x€B \ |xo + ei| / 

since 2^^^^ — Ci G SM^. Hence the function v - extended trivially on 9B \ {— ci} - is 

continuous in i? \ {— Ci}. We fix a unit vector e G perpendicular to Ci (i.e., |e| = 1 and 
e • Ci = 0), and we show that v is symmetric with respect to the hyperplane T := {x G M'^ : 
X ■ e = 0}. For this we apply a moving plane argument based on the integral representation 

(4.3) and reflection inequalities derived in [4,11] for Gi. We need some notation. For A > 0, 
we consider the open half-space Hx = {x G : x ■ e > A} and the reflection x i— > x'*' := 
X — 2(x ■ e — A)e at the hyperplane dHx. We also consider the set 

Ja := {x G B : x ■ e < A and x'^ ^ B} 

which has nonempty interior if A > 0. With these definitions, the inequalities stated in [4, 
Lemma 4] (see also [11, Lemma 3] for the biharmonic case) translate into the following 
reflection inequalities: 

Gi{x\y^)>Gi{x,y^) and 

(4.4) , , , , } for all X, y G -f^A n B 
G^{x\y')-G,{x,y)>G^{x,y')-G^ix\y) ^ 
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and 

(4.5) G'i(x\y) - > ioTxeHxnB,ye -h. 
Now (4.5) and the strict positivity of f in B imply that 

(4.6) / [Gi{x^,y) -Gi{x,y)\h{y,v{y))dy> Q for A > and x G i^A n B. 

We claim that the following reflection inequality holds for every A > 0: 
(Ca) v{x) < v{x^) for all x E HxHB. 

We put 

A* := inf{A > : (Cv) holds for all A' > A}. 

Then A* < 1. Using the continuity of v in B, it is easy to see that (CaJ holds. We suppose 
for contradiction that A,, > 0. Since < |x + ei|~° < \x^* + ei|~° for x G Hx^ fl B and v is 
positive in B, (CaJ yields 

(4.7) h{x,v{x)) < h{x^%v{x^*)) for x G i^A, n B. 
We claim that 

(4.8) v{x) < v^x^*) for all a; G i^A, n B. 
Indeed, by using (4.4), (4.6) and (4.7) we have 

v{x^*)-v{x)= / [Gi{x^%y) - Gi{x,y)]h{y,v{y))dy = / . . . dy + . . . dy 

= [ ([Gi(x^* , y) - Gi(x, y)] h{y, v{y)) + [Gi(x^* , /*) - Gi(x, y^')] h{y^%v{y^'))) dy 
+ / \G^{x^\y) - G^{x,y)\h{y,v{y))dy 



> / [Gr{x'\y'*)-G^{xy')][h{y'^%v{y'*))-h{yXy))]dy>Q for x G i^A, n B. 

J Hx^riB 

Hence (4.8) is true. For < /i < A* we now consider the difference function 

w^:if^nB^M, w^{x) = v{x^) - v{x) 

and the set 

W^:= {xeH^r^B ■ w^{x) < 0}. 
We note that = 0, and we claim that 
(4.9) \Wf,\ ^0 as /i ^ A„ 

where | ■ | denotes the Lebesgue measure. Indeed, let e G (0, A^,), and consider the compact 
set 

K := {x eB : x ■ e > A* + e, |x| < 1 - 5} C i^A. H B 
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Then inf Wx^{x) > by (4.8). The continuity of in B imphes that there exists Ai G 
(A* — e, A*) such that 

inf wJx) > for Ai < /i < A*. 
Hence C {H^ n B) \ ii' C {x G B : |xi — A*| < e or |a;| > 1 — e} and therefore 

\Wfj,\ < 2^e + euN-i for \i < fi < K, 

where wn^i denotes the area of the — 1-dimensional unit sphere. Since e was chosen 
arbitrarily small, (4.9) follows. Next we note that, for fi > ^ and x G W^, 

hix^'^vix^')) -h{x,v{x)) = \x^ + ei\-''v''{x^) - |x + ei|~"t;'?(x) > |x + ei|"°(t;«(x'^) - ^;«(x)) 



(4.10) 
with c(A* 
(4.11) 



> 



A. 



'[v'^ix") - v'^ix)] > (y) qv''~\x)[v{x'') - v{x)] > c{X,)w^{x) 



q[ sup v[x'^ 



9-1 



< oo. We also note that 



h{x^, vix^")) - h{x, v{x)) > ioix^H^X W^. 

From now on we assume that ^ < < A*. For x G we use (4.4), (4.6), (4.10) and (4.11) 
to estimate 

> w^{x) = 



> 



> 



[Gi(x^, y) - Gi(x, y)] h{y, v{y)) dy 
[Gi(x^ y) - Gi(x, y)] h{y, v{y)) + [Gi(x^ y^) - G,{x, y^] h{y^, v{y^))]dy 
[Gi(x^ y^) - Gi(x, y^)] [/^(|/^ v{y^)) - h{y, v{y))] dy 
> [ [G,ix^,y'^)-G^ix,yn][hiy'',viy'^))-hiy,viy))]dy 

JWa 



>c(A,) / [G^{x'',y'')-Gi{x,yn]w,iy)dy 

(4.12) >CN,mc{K)j \x-y\^~^w^{y)dy, 

where in the last step we use the estimate 

< Gi{x, y) < CN,m\x - y\^~^ ioi x,y eB, x ^ y 

with some C]^^rn > 0, which is easily deduced from the integral representation of Gi in 
Section 3. Next we pick s > 1 large enough such that 1 < q := xzx ^ ^- Then the 
Hardy-Littlewood-Sobolev inequality (see e.g. [17, Section 4.3]) implies that 



(/ 


/ \x - y\^ ^Wf,{y)dy 


S \ s 

dx 1 = 


1 ■ T"^ * (w.t^M) 











L" 



q\\W^,\\L^i{W^,) 
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with a constant Cs^q > 0. Combining this inequahty with (4.12) and Holder's inequahty, we 
obtain 

(4.13) ||w^f.||L»{iyp) < Co||u;^||l9{w/^) < Co|W^|~||ui^||L^(iy^) with Cq := C7v,„, c(A*) c^,,. 

Now (4.9) and (4.13) imply that = if yU < A* is close enough to A*. Hence 

property (C^) holds if /i < A* is close enough to A*, which contradicts the definition of A*. It 
follows that A* = 0, thus (Cx) holds for all A > 0, as claimed. By continuity, we now deduce 
that 

v{x) < v{x°) for X e HoHB. 
Repeating the moving plane procedure for — e in place of e, we get 

v{x) > v{x°) for X e HoHB. 

Hence equality holds, i.e., v is symmetric with respect to the hyperplane T = {x G : 
X • e = 0} as claimed. Since e was chosen arbitrarily with |e| = 1 and e ■ ei = we conclude 
that V is axially symmetric with respect to the xi-axis. □ 

Corollary 13. The function u only depends on the xi-variable. 

Proof. Since v is axially symmetric with respect to the Xi-axis, the same is true for u. Let 
z G M^~^ be arbitrary, and consider the function 

U,:R^^ M, U,{xi, x') = u{xi, x' - z) for (xi, x') G [0, oo) x R^~^. 

Then Uz satisfies the same assumptions as m, so it is also axially symmetric with respect to 
the xi-axis. This readily implies that u only depends on the xi-variable. □ 

Theorem 14. Let f : [0, oo) — > [0, oo) he a continuous function with /(O) = and f{s) > 
for s > 0. If u is a classical non-negative bounded solution of the one- dimensional problem 

(_l)-^(2'n) = j(^) (0, oo), m(0) = u'iO) = ...= u^"'-'\0) = 
then u = 0. 

Proof. The differential equation admits a first integral given by 

m—1 ^ 
H = ^(-l)*M«M(2m-*) ^ (_i)mj^_(^(m)^2 ^ ^(^^^^ 

i=l 

where F{s) := Jq f{s) ds. Indeed, we calculate that 

^ = ^(-1)^(^^(^+1)^(2—0 +^W^(2m-(i-l))j ^ (-1)-(^m("^)m(™+^) + f{u)u'^ 
1=1 

= -n'n(2"^) + (-1)— iu('")u(™+^) + (-l)"^(^u('")u('"+^) + f{u)u'^ 

= n'(^(-l)'"/(n) -^(2"*)) = in(0,oo). 

Suppose for contradiction that u ^ 0. Since u has a Green function representation by 
Theorem 9, we infer that u is strictly positive in (0, oo), so that ^(2™) = (— 1)('")/(m) has no 
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zero in (0, oo). By the mean value theorem, this imphes that u^^^ has at most 2m— j zeros in 
(0, oo) for j = 0, . . . , 2m. Hence every u^^^ is eventually monotone and has a limit as t — >■ oo 
since it is bounded. From this it clearly follows that 

u^^\t) ^0 as t ^ oo for j = 1, . . . , 2m, 

but then also u[t) ^ by using the equation and the assumptions on / again. Since -F(O) = 0, 
we thus find that i7 = in [0, oo). In particular, the boundary conditions yield 

o = i^(o) = ^^[«(-)(o)]^ 

so that n^™''(0) = 0. If m = 1 we conclude m(0) = u\Q) = 0, so that n = by the uniqueness 
of the solution of the initial value problem - contrary to what we have assumed. 
If m > 1, we set v = —u" and find that is a bounded solution of 

(_l)™-i^(2(— D) = > in (0, oo), v{<S) = t/(0) = . . . = ?;('"-2)(0) = 0. 

So V also has a Green function representation by Theorem 9, and thus f > in (0, oo). In 
sum, we have 

m(0) = and m > 0, m" < in (0, oo), 

which forces m'(0) > and therefore contradicts the boundary conditions. The proof is 
finished. □ 

Proof of Theorem 4 (completed). By Corollary 13, we have u{x) = u{xi) with some 
function u : [0, oo) — > M. Since u satisfies the assumptions of Theorem 14, we conclude that 
u = and therefore u = 0. □ 

5. Appendix 

Lemma 15. Let v be a classical bounded solution of {—A)"^v = g{v) in M^. If g : M — [0, oo) 
is convex and non-negative with g{s) > for s < then v > 0. 

Proof. Notice that the boundedness of v implies the boundedness of D^v for j = 1, . . . , 2m. 
First we show that (— A)'i; > in for / = 1, ... ,m — 1. Assume that there exists / G 
{1,. . .,m-l} andxo G with (-A)'w(a;o) < but {-AYv > in for j = / + 1, . . . , m. 
We may assume w.l.o.g. that xq = 0. Let vi := (— A)'w for Z = 1, . . . , m — 1 and set vo = v. 
Then we have 

-Avo = vi, -Avi = V2, ... - Avm-i = g{vo) inR^ . 

If we define spherical averages w{x) = ^iv-\^^ (q) w{y) day, r = \x\ then the radial func- 
tions Vq, vi, . . . , Vm-i satisfy 

-Avq = vi, -Avi =V2, ... - Avra^i > g{vo) in R^, 

where we have used Jensen's inequality and the convexity of g. Since vi{0) < we also have 
vi{0) < 0. Moreover {r'^~^v'i)' = —r'^~^vi^i or < —r'^~^g{vo), which in both cases is non- 
positive. Since v'i{0) = we see that v'i{r) < for all r > 0. In particular vi{r) < vi{0) < 0. 
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Next we integrate the equation 

Avi^i = -vi > -vi{0) > 

and obtain r^~^v'i_j^{r) > —^vi{0), i.e, vl_i{r) > —j^vi{0). The unboundedness of v'i_i 
yields a contradiction. 

Finally, we need to show that v = vq > 0. Assume that ^0(2^0) < and w.l.o.g. Xq = 0. 
Since Avq = —vi < we see that v^ir) < and define a := linir^oo ^o('") < 0. Then 
limr^oo Av„^_i{r) = —g{a) < 0, i.e., Avm-i{r) < —■^g{a) < for r > r^. By integration this 
leads to linir^oo = —00, which contradicts the boundedness of w^-i- This finishes 

the proof of the lemma. □ 

Acknowledgement. The authors thank Hans-Christoph Grunau (Magdeburg) for helpful 
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